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ABSTRACT 

We investigate coronal heating and acceleration of the high- and low-speed solar wind 
in the open field region by dissipation of fast and slow magnetohydrodynamical (MHD) 
waves through MHD shocks. Linearly polarized Alfven (fast MHD) waves and acoustic 
(slow MHD) waves travelling upwardly along with a magnetic field line eventually form 
fast switch-on shock trains and hydrodynamical shock trains (N-waves) respectively 
to heat and accelerate the plasma. We determine one dimensional structure of the 
corona from the bottom of the transition region (TR) to lAU under the steady-state 
condition by solving evolutionary equations for the shock amplitudes simultaneously 
with the momentum and proton/electron energy equations. Our model reproduces the 
overall trend of the high-speed wind from the polar holes and the low-speed wind from 
the mid- to low-latitude streamer except the observed hot corona in the streamer. The 
heating from the slow waves is effective in the low corona to increase the density there, 
and plays an important role in the formation of the dense low-speed wind. On the other 
hand, the fast waves can carry a sizable energy to the upper level to heat the outer 
corona and accelerate the high-speed wind effectively. We also study dependency on 
field strength, Bq, at the bottom of the TR and non-radial expansion of a flow tube, 

/max, to find that large -Bo//max ~ 2 but small ~ 2G are favorable for the 
high-speed wind and that small B^/ f-aia^ ~ 1 is required for the low-speed wind. 
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1 INTRODUCTION 

The problems of the coronal heating and the solar wind ac- 
celeration should be investigated simultaneously since they 
are very closely related, although they have been often dis- 
cussed separately. Even combined models, taking into ac- 
count them together, usually adopt too simplified heating 
and acceleration law: they adopt unspecified mechanical 
energy and momentum inputs of the exponential type, oc 

(r-RQ), 



exp( ; 

lin, for lack of an alternative (e.g.Sandbak 



, with an assumed constant dissipation length, 

Leer 1994), al- 



though the constant is poorly supported by fundamental 
physical processes. Therefore, it is required to treat the en- 
ergy transfer and momentum transfer simultaneously in the 
broad corona including the solar wind with realistic heating 
and acceleration functions. 

The origin of the energy to heat the corona and 
accelerate the solar wind is generally believed to lie in the 
turbulent convective motions beneath the photosphere. The 
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major difficulties in understanding the problems are how to 
hft up the energy to the corona and let it dissipate at the 
appropriate locations. Various modes of waves are supposed 
to play important roles in these processes (see Roberts 2000 
for recent review) . Not only the surface granulations but also 
dynamical a ctivities such as the small magnetic reconnec- 
tion events ('Tsuneta"l996!; 'TarbeU, Rvutova, & Covingtonl 
[1999; Sakai ct al., 2000,; Katsukawa & Tsuncta 2001) ex- 
cites th e waves from a t the photospheric level to in the 
corona iSturroclJ 119991) . MHD waves have been studies 
long as reliable heating and acceleration sources of the 
solar coronal plasma (e.g. Barnes 1969; Belcher 1971). 
Many groups have examined advanced properties on the 
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waves, and shock formation ( H.oUweg 1992) to study the 



wave heating. As an effect bey ond the MHD, resonant dissi- 
pation of ion-cyclotron waves jDusenberv fc Hollwedll98ll : 
iMarsch. Coertz. fc Richtedll98a) " has also been highlighted 
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as a h eatinK source of the heavy ions i n the hiRh-speed sola r 
wind dCranmer. Field, fc Koh]ill999l : IHu fc HabbaJ [l999l) . 
However, because of the too high resonant frequencies, the 
protons and electrons, underlying the coronal plasma, are 
expected to be mainly heated by the MHD waves with lower 
frequency in which the most of the wave energy is con- 
tained. Signatures of such MHD waves have been observed 
in various port io ns of the solar coron a ( Hasslcr ct al. 19901 
|BarMne^£^^J_ Il99^: lOfman Nakariakov. fc DcD orestI 
I1999I: lO'Shea. Muldach. fc Fleck! 120021) . whereas dominant 
processes in the wave heating yet remain to be determined. 

When one considers the MHD waves travelling along 
with a magnetic field line in the low-/3 plasma, they can be 
divided in two types in the linear regime where the wave am- 
plitude is sufficiently smaller than the phase speed. One is 
slow-mode wave which is longitudinal wave and essentially 
identical to the acoustic wave. The other is Alfven wave 
which is transverse wave. When the non-linear terms are 
taken into account, the degeneracy of the transverse waves 
is removed; they can be divided into circularly polarized 
Alfven wave classified as the MHD intermediate-mode, and 
linearly polarized Alfven wave classified as the MHD fast- 
mode. Among these three types of the MHD waves, the 
fast-mode and slow-mode waves nonlinearly steepen their 
shapes through propagation to eventually form the MHD 
fast shocks and slow shocks respectively. Suzuki (2002; here- 
after S02) investigated the dissipation of the N-wave which 
is a special case of the slow-shock trains propagating along 
the field line. S02 also inspected the consequent heating in 
self-consistent modelling to conclude that the N-waves can 
heat the inner corona effectively though they cannot heat the 
entire corona. The heating by the dissipation of the switch- 
on shock trains, a special case of the fast shock trains, was 
proposed by HoUweg (1982; hereafter H82). He studied the 
problem on fixed coronal structures to conclude that it can 
be a reliable process of the coronal heating. However, the 
self-consistent treatment of their propagation and dissipa- 
tion has not been performed so far. 

In this paper, we investigate the coronal heating and the 
solar wind acceleration by the switch-on shock trains as well 
as the N-waves in a self-consistent manner. First, by extend- 
ing the formulation introduced by H82, we derive an equa- 
tion describing variation of the switch-on shock trains in the 
similar way to that adopted bv lStein fc Schwart j il972l) and 
S02 for the N-waves. We solve these equations for the shock 
trains simultaneously with the equations of the momentum 
transfer and energy transfer to determine the unique struc- 
ture of the corona and solar wind from the bottom of the 
transition region (TR) to lAU for given wave energy fiux 
of the linearly polarized Alfven waves, i^w,±,o, and acoustic 
waves, F„j|_o, at the coronal base. The heating and the ac- 
celeration by the dissipation of the shock trains are explicitly 
taken into account. Therefore, we do not have to take the 
phenomenological approach adopting the exponential type 
of the heating law. Conversely, we can test the validity of 
the assumption of the constant dissipation length for our 
processes. 

We further explore possibilities whether our model can 
account for the difference between the high-speed solar wind 
( ~ 700km/s) from the polar coronal holes and the low- 
speed solar wind ( ~ 400km/s) in the equatorial region 



during the low solar activity phase. We adopt various field 
strength at the bottom of the TR and non-radial expansion 
factor as well as Fw,±,o and i^w,||,o to study which parameter 
plays a role in exhibiting the distinctive difference between 
the two types of the solar wind. 



2 VARIATION OF SHOCK AMPLITUDE 
2.1 N-waves 

Acoustic waves have not been considered as a major heat- 
ing source of the solar corona for several decades, be- 
cause those generated by the granule motions at photo- 
spher ic level cannot reach the corona iStein fc SchwarO 
119721) . However, recent observations have revealed dynami- 
cal natures of the corona; a lot of transient events such as 
microfiares constantly occur. Some models JSturrocklll99gl : 
iKudoh fc Shibata 1999) describing those events predict ex- 
citation of acoustic waves, or equivalently slow-mode MHD 
waves, not at the photosphere but in the corona. Moreover, 
Ofman et al.(1999) identified slow MHD waves propagating 
in polar plumes. If one considers such waves travelling along 
with the magnetic field line in the open flux tube region, 
they eventually form shocks and propagate upwardly as the 
N-waves (Fig0, a special case of the slow-shock trains, with 
heating the surrounding plasma. For example, the waves of 
sinusoidal shape with period, ry ~ 100s, steepen their wave 
fronts to form the shocks after propagation by a distance 
of ~ lO^km in the corona, provided that the initial wave 
amplitude is ~ 10% of the sound speed (S02). 

A model for the N-wave s was developped for pro paga- 
tion in the chromosphere bv lStein fc Schwartd il972h . and 
generalized for propagation in the corona by S02. Now, we 
briefly summarize the derivation of an equation for variation 
of amplitude of the N-waves by following S02. We assume 
the shock amplitude, ^Wshj, is sufficiently smaller than am- 
bient sound speed, Cs (weak-shock approximation). Then, 
entropy generation at each shock front is given by the jump 
condition across the shock as 



Aaii _ 7-H Ap 3_ 27(7 -H) 3 
R 1272 ^ p ' 3 * 



(1) 



iLandau fc Lifshitdll959l) . where R is the gas constant, 7 
is a ratio of the specific heat, Ap denotes pressure differ- 
ence at the front, p is pressure in the upstream region, and 
Qgh,||(= 5ah,|i/cs) is the normalized shock amplitude. Using 
the perfect gas law, T = p/{pR), the energy loss rate —Q\\ 
per volume (erg cm~^s~^) at the shocks for the waves with 
period ry is derived from eq. 



— — poTAs||/r|| — 



27(7+ l)Poai.|| 



3rii 



(2) 



Divergence of wave energy, E\^^ = ■^■ypal^^ ||A||, per wave- 
length, A||, for a shape constant, <;||(= 3 for the N-waves 
illustrated in Fig0, can be written as 



(3) 



where A is a cross section of a flow tube. Substitution of 
eq.J^ into the logarithmic derivative of Ex^^ gives variation 
of the shock amplitude in the static media as follows (S02) : 
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Figure 1. Schematic picture of tlie steepening of acoustic waves 
(top) and linearly polarized Alfven waves (bottom). 



da. 



dr 



Qsh.ii , 1 dp -^ii 2(7 + l)ash, 
2 p dr 3 CsTi 



1 dA 3 dcs 
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This equation determines the variation of the shock ampli- 
tude of the N-waves in the solar corona according to the 
physical properties of the background coronal plasma. The 
first term in the right-hand side denotes the amplification in 
the stratified atmosphere, and the second term indicates the 
loss by the shock heating. In the (1 — 2) J?© region, where the 
dissipation is important to heat the corona, these two terms 
always dominate the other terms, the third term arising from 
the geometrical expansion and the fourth from temperature 
variation. 



2.2 Switch-on Shock Trains 

The various dynamical activities on the solar surface are 
expected to continuously produce transverse waves. Some of 
them would be linearly polarized Alfven waves propagating 
upwardly. Phase speed of the linearly polarized Alfven waves 
varies as a function of fluctuatin g field, SB±, normal to the 
background field, _B|| iKulikovskiv &: Lvubimov^ll965^ : 



(5) 



The above equation al so satisfies the req uirement of the 
plasma compressibility iMonteomervlllQsil) . 

{Bl + 6Bl)^^^ 



const. 



(6) 



Equations ^ and @ shows that both the wave crest and 
the wave trough overtake the wave front with a speed. 



52 + SB^, 



Bn 



VA[ 



Bf^+SBl 



1], 



(7) 



where we here define Alfven velocity, da = B^^/y'^npo for 
the ambient density po at the wave fronts. Then two shock 
fronts are formed per wavelength as illustrated in Fig. Q 
which is also seen in results of time-dependent calculations 
by Nakariakov et al.(2000). The figure also displays that as 



moving from the upstream region to the downstream re- 
gion, 5B± and 5v± 'switch-on' from zero to finite values. 
Therefore, it is called switch-on shock, a special case of the 
fast shocks. Periodic trains of the switch-on shocks prop- 
agate upwardly with a speed ~ va, similarly to the usual 
Alfven waves except that they dissipate the wave energy at 
each shock. 

Let us estimate a distance the sinusoidal wave travels 
till the formation of the shock fronts. We can categorize 
these transverse waves into two types in terms of location of 
the wave generation. First type is the wave generated by the 
granule motions at the photosphere. Since the Alfven veloc- 
ity in the photosphere and low chromosphere is quite small, 
the sinusoidal waves steepen to form the shock trains before 
reaching the corona. This is seen in re cent time-dependent 
calculation bv lKudoh fc Shibatal 1^9^, illustrating that the 
transverse waves created at the photosphere propagate into 
the corona as the fast shocks. 

Second type is the wave excited by transient events such 
as microfiares in the corona. To form the front, the wave 
crest or trough have to overtake the wave front by a quarter 
of a wave length, Aj_ . If we set that the sinusoidal waves are 
generated at time to and position Rq, time t and position R, 
at which the waves steepen to form the fronts, satisfy the 
conditions below in a static atmosphere. 



dt va[(- 



B^. 



-1] 



dr[{- 

Ro 



B^ 



If we consider wave with period, t± = 120s, generated at 
-Ro ~ Rq -I- 1 X 10'*km, substitution of typical values in the 
corona of ~ 5G, SB± ~ IG, and density scale height, 
H ~ 6 X 10''(j^)km, into eq.® gives R ~ I.SRq. It in- 
dicates that even if the linearly polarized Alfven waves are 
excited in the corona, they form the fast shocks in the low 
corona. Moreover, the assumption of the generation of the 
sinusoidal waves is idealistic and the initial wave shape is 
supposed to be more or less out of shape in the real solar 
corona. In a such case the distance till the formation of the 
shocks should become shorter. 

From now, we would like to derive an equation describ- 
ing variation of amplitude of the switch-on shock trains by 
extending the formulation by H82. From the jump conditions 
across the switch-on shock front (Boyd & Sanderson 196iD, 
entropy generation As± at the shock can be expressed as 



8np 



If 



(9) 



for the weak-shock condition of cr — 1 ^ 1 (H82), where p is 
pressure in the upstream region and cr is a ratio of the down- 
stream density to the upstream density. Field strength, SB^, 
in the downstream region perpendicular to the underlying 
field, Bii , can also be derived from the shock conditions: 



5Bi 



2(c 



(10) 



where Cs and va is the sound speed and Alfven speed in the 
upstream region. Eliminating (a — 1) from eqs.Q and llOII . 
As_L is rewritten as a function of 5B± or 5v±: 

B^ 



327rp(l 



327rpi(l-ci/i;i)2' va 



.^B±^4 
B|| ' 

^ 5v±_ ^4 



(11) 
(12) 
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r(ash,j 



(13) 



where we have used the relation of Sv± = Sb± / ^A-Kp in the 
downstream region. Recalling that the two shock fronts form 
per wavelength in this case (H82; Fig0, energy loss rate 
—Q±{erg cm~"^s~^) for the waves with period t± becomes 



1 4 

= -PT^S± (-T^) = -— ^ 

2 167rrx (1 - c|/i;^)^ 



(14) 



It should be noted that the heating, Q±, is proportional 
to Qgij X , which is a higher order term with respect to cvsh 
compared to the case of the N- waves giving oc af^^ y 
(eq.l|5J), since the liberated energy at the switch-on shocks 
is transfered to magnetic field as well as heat. This fact 
implies that the damping by the switch-on shocks is weaker 
and the periodic trains could travel to a more distant region 
than the N-waves to contribute to the heating in the outer 
region. 

Now we derive the equation for Qsh.x by using eg. 11411 
under the WKB approximation. In moving media with ve- 
locity V, even though the dissipation does not occur, the 
wave energy flux is not conserved. Instead, one can have 
a new cons erved quantity of wave action constant which is 
defined as (lJacaues|ll977ft 



pa'^i^j_VA(v + va){v + va) 



(15) 



where ?x is a shape constant for the switch-on shock trains. 
In our calculations we set ^± = 3, corresponding to the wave 
illustrated in Fig0 Although Sw,_l is the conserved quan- 
tity without dissipation, we have to take into account the 
energy loss at the shocks for the switch-on shock trains. 
Then, an equation governing variation of Sw,_l becomes 



VA 



-V ■ Sw.i = 



B 



VA+v levrrx (1 - ci/i;i)2 ■ ^^^^ 

Since the left-hand side of eg. (1161 is rewritten in terms 
of wave energy flux, i^w,_L, and wave pressure, Pw,x, as 



jV • S-^,± = V ■ Fw,x — V ■ Vpw.x, eg. ljl^ll essentially 
indicates variation of the wave energy flux. In eg. (1161 wave 
period, rx, is measured from the observer moving with the 
flow. It is more useful to consider the period, 



t±va/{va + y), 



(17) 



as seen by the stationary inertial observer, since it is a con- 
stant along with a given field line (H82). Therefore, eg. 1161 
is rewritten as 



Bi 



IG^rx.i (1 - cl/vi)^ ' 



(18) 



where we here assume S^,± has the only radial variation: 



1 d(^S„ I ) 



S'w .X — \S\f 



derivative of S, 
1 dS, 



dr 



llSl ) gives 
2 dosh, 



The logarithmic 



p dr Qsh.x dr 

3vA + y dVA _|_ 
VAivA+v) dr VA-\-vdr 



dv 



(19) 



Then, an equation describing variation of the shock ampli- 
tude in the steady fiow is derived from eqs. 1181 and 1191 : 



das 



Qsh, 



dr 



1 dp 
p dr 



?xB||Qsh._ 



8rx,iv^(uA + u)2(l - ci/viy 



1 dA 3va + V dvA 
A dr va{va + f) dr 



dv. 



Va + V dr 



(20) 



We would like to emphasize that eg. 1201 is derived for the 
first time, whereas eg. 1181 is directly solved in the fixed 
background structures in H82. An advantage of using eg. 1201 
instead of eg. 1181 is that we can directly get the shock 
strength, ash.x, when performing the numerical integration. 

In the static limit {v <^ va), eg. 1201 shows an analogous 
form to eq.lQJ for the N-waves: 



da. 



'ah,. 



dr 



1 dp 

p dr 4:T±va{1 - 



-ci/vir 



1 dA 3 dvA , , . 
A dr VA dr 



In the square brackets of the right-hand side, the 1st and 
3rd terms are identical to those in eq.l0J, and in the 4th 
term va appears instead of c^. An important difference is 
seen in the 2nd term denoting the dissipation at the shocks. 
The variation of the amplitude is proportional to a^j^ in 



the case of the switch-on shocks, while that to a„ 



in the 



N-wave case, which reflects the different dependencies of the 
entropy generation on ash (eqs.Q & 1131 ). This also shows 
that the switch-on shock trains are less dissipative. 



2.3 Limitations in Our Modelling 

We had better comment on the limitation when using eqs.Q 
and 1201 for the analysis of the propagation of the shock 
trains. An interaction between the N-waves and the switch- 
on shock trains is supposed to give little influence, since it 
is purely the non-linear effect. Our calculations show that 
the shock strength is flnite but still small (ash < 1) (©• 
Consequently, the two modes of the shock trains can travel 
almost independently, hence, it is justified to solve eqs.Q 
and lUnj 

separately. 

The first limitation on the propagation of the N-waves 
is that it is valid only in the static media. However, this 
is accomplished in our calculation since the N-waves with 
r|| ~ 100s mostly dissipate in the low corona where the 
static approximation is sufficient (33 ■ The second limita- 
tion is that it does not take into account the effect of the 
gravity. This is also negligible for the waves with ~ 100s, 
because the acoustic cut-off period for the acoustic gravity 
wave is much greater (rac ~ 1000s) in the corona. The third 
limitation is that ash,|| should be less than 1, since eq.(|lj 
is derived under the weak shock approximation. As shown 
in the results, our model does not break this limitation ei- 
ther, whereas the condition at the coronal base is almost 
marginal, Qsh,|| — 0.8 (FigEJ. The final limitation is that we 
do not take into account viscosity explicitly. Although we 
consider viscosity implicitly through the shock, compressive 
viscosity, which is not negligible in the corona, may make 
the acoustic waves dissipate even when they mildly steepen 
before the shock formation (Ofman et al.2000). However, 
this effect dose not affect the coronal structure, since the 
dissipated energy from the acoustic waves is mainly lost as 
the downward thermal conduction, unrelated to how rapidly 
they dissipate (see 

As for the switch-on shock trains, we do not have to 
consider the problem in the static media as eg. 1201 has been 



Coronal Heating and Solar Wind Acceleration by MHD shock trains 5 



derived with taking into account the steady flow. The effect 
of the gravity is negligible, since the switch-on shock trains 
are basically transverse. The condition of the weak shock, 
Qsh,± < 1, is also accomplished (FigBJ- The necessary con- 
dition for the WKB approximation is that the scale height of 
the Alfven speed. Ha, is sufficiently larger than the wave- 
length. The condition is easiest to be violated in the low 
corona where Ha ~ 2 x lO^km and va ~ 10"^km/s. There- 
fore, the WKB approximation becomes marginal there for 
the waves with t± ~ 200s, while it is reasonable for those 
with T±_ ~ 100s. In more sophisticated models dynamical 
treatments (e.g. Ofman & Davila 1997; 1998) are required 
for those low-frequency waves. The most severe limitation 
is that we assume that the transverse waves dissipate only 
through the switch-on shocks. Although taking into account 
multi-dimensional effects is beyond the scope of this pa- 
per, transverse variation of the f ield strength leads to the 
dissipation by the phas e mixing jHevvaerts fc Priestlll983 
ISakurai fc Graniklll984ri and through the mode conversion 
ij^ike^eto^^ob^s^!^ Murawskilll99'^ . If the waves are 
not coherent and show more or less turbulent-like natures, 
the turbulent cascade might occur (e.g. HoUweg 1986; Hu, 
Habbal, & Li 1999). To determine the dominant process(es) 
in the solar corona and the solar wind is an important task. 
By comparing our results with those considering the other 
processes, we can proceed a further step to understand the 
coronal heating and the solar wind acceleration. 



3 MODEL FOR CORONA AND SOLAR WIND 
3.1 Basic Equations 

In this study, the solar corona and wind are assumed to be 
fully ionized MHD plasma consisting of protons and elec- 
trons. The fluid (MHD) treatment might be controversial 
especially in the outer region (r ~ IO-R0) as the plasma 
becomes coUisionless with respect to the Coulomb collision, 
and the kinetic treatment would be required in more realis- 
tic models. However, randomness of the magnetic flelds may 
operate to let them exhibit the MHD characteristics. 

We here present basic equations to describe coronal 
wind structure in a flow tube with a cross-section of A under 
a steady state condition. The equation of mass conservation 
is written as 



pvA — const., 



(22) 



where mass density, p, is related to proton and electron 
number density, rip and rie, and proton and electron mass, 
rrip and rric, as p — npirip + 12^171^ — ripTTip. We assume 
Tic = Up = n, and the proton and electron components have 
the equal outflow velocity in the solar wind. Then, the equa- 
tion of momentum conservation becomes 



dv 
dr 



GMp 



Idp 

p dr 



dr 



dr 



(23) 



where Pw.x is the wave pressure of the switch-on shock trains 
and p„ j| is that of the N- waves. Gas pressure, p, is related 
to p and temperature of the protons, Tp, and the electrons, 
Te by an equation of state for an ideal gas; 



P ■ 



-{Tp + n: 



(24) 



where fce is the Boltzmann constant. 

As for the energy transfer, the protons and electrons are 
allowed to have different temperature. The energy equation 
for the protons is 



n"' d .k^Tp, 



7 — 1 dr n^-i 
and that for the electrons is 



= -V ■ Fc,p + Cpe -I- Qx,p + I 



n? d , ksTc 

1> — ( T 

7-I dr n-i-^ 



-Ccp+Q±,c+t 



(25) 



2$(Te),(26) 



where Cpc ~ —Cop indicate energy exchan ge between 
the electrons and protons by Coulomb collision iBraginskii I 
1965; Sandback & Leer 1994), Q±^p and Qx,c denote heating 
from the switch-on shock trains to the protons and electrons, 
Qll^p and Q\\^e axe heating from the N- waves to the pro- 
tons and electrons, and $(Te) is r adiative loss function by 
iLandini fc Monsignori-Fossi lll99(Jl . fc,p and Fc,e is thermal 
conduction taken by the protons and electrons, whereas we 
here use the classical Spitzer conductivity by the Coulomb 
collisions: 



F ^ 

^ c.p — f^p J 

dr 



5 dT 
Kp,0Jp , 



dTc 
dr 



dr 



(27) 



(28) 



where n^fi — 3.2 x 10~^ and He,o = 7.8 x 10~^(erg 
cm~^s"^K"'^''^) jBraeinskii Ill965^ . It should be noted that 
adoption of the above classical forms of thermal conduction 
is controversial especially in the outer corona and the so- 
lar wind region, since the plasma becomes coUisionless there 
(e.g. HoUweg 1976). 

In the very inner corona where the coUisional coupling 
between the protons and electrons is satisfied, we use one- 
fluid approximation for the energy transfer. The energy 
equation for the one-component plasma can be obtained by 
combining eas. l25t and I26II : 



2nT d , ksT , 

V — ( 

7-I dr nT-i' 



-v-Fc + Qx + 0||-n ^(rc: 



(29) 



where T = {Tp + T,)/2, ~ -(Kp,o + ^3.0)^2^:, = 
Qx,p + Qx,o, and Qy = Q\\^p + Q\\,c- Please note that total 
number density becomes 2n(= rip -I- rio). A practical way to 
connect the coronal wind structures of the two-component 
plasma to those of the one-component plasma will be de- 
scribed in next subsection. 

For the wave pressure terms in ea. (l23ll . we take 
the WKB approximation. Wave pressure of the switch-on 
shock trains is essentially identical to that of the usual 
Alfven waves. Then, the term as to the wave pressure gra- 
dient in eg. 1231 can be written as 



dp. 



dr dr 87r?x 



^ d QsVx-By , 
dr ^ 87r?x ' 



(30) 



Although the wave pressure term for the N- 
waves can be found in S02, we neglect it since inclusion of 
this term does not affect the results at all (S02): 



dr 



= 0. 



(31) 



The heating terms, Q± and Qy, have been already derived 
in egs. 11811 and Q. Distribution of the dissipated energy at 
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the shocks to the protons and electrons is unclear. Therefore, 
we use following parameterization : 



fex/ii)Q±/ii 



Q 



ix/ll^i-i/ll 



(32) 
(33) 



We perform our calculation for two extreme cases of /ix/|| = 
0.5 and 0, whereas we mainly discuss model results adopting 



/II 



0.5. 



Since the switch-on shock trains travel in the very sim- 
ilar manner to the usual Alfven waves except the shock dis- 
sipatio n, wave energy flux, i^w,x, is explicitly expressed, fol- 
lowing |jacaue3 1^3), as 



-pash,xt'A[i'A(l + ash,±)* + o'"] 

— pash,x«A[uA -I- \v]. 
?x ^ 



(34) 



For the N- waves, we use a form in the static media (S02), 
-Fw.ll ^ — pa^h.iiCsCs(l + i ^ «Bh,||), (35) 



1 2 2 1,^1 + ^ 



since they dissipate quickly ( H5.2II and work effectively only 
in the low corona where the static approximation is fulfilled. 

To take into account the non-radial expansion of the 
flow tube due to conflgurations of the magnetic field, the 
cross-sectional area, A., is modeled as 



A = r 



where 



2 /"maxG 



(i — ri)/<7 



+ .fl 



(36) 



/l = l-(/n.ax-l)e(l-'-^'/^ 

JKopp fc Orralll976l:IWithbroell98ll . The cross section ex- 
pands from unity to /max most drastically between r = ri — a 
and ri -I- a. Of the three input parameters, /max is the most 
important in determining the solar wind structure. In this 
paper, we consider cases between /max ~ 1 and /max ~ 20. 
As for the othe r two parameter s, we employ ri = 1.25ii0 
and a = O.IRq JWithbrodllQSll . 

The geometrical expansion of the flow tubes is subject 
to the open magnetic field lines. Conservation of magnetic 
fields, V ■ i? = 0, gives a condition for radial magnetic field, 
B\\ , as 



B\\A = const. 



(37) 



In our modelling, we leave the radial field strength, Bq, 
at the inner boundary located at the bottom of the TR 
as an input parameter, and determine -B|| (r) according to 
eg. 13711 . Note that Bo is supposed to be smaller than the 
field strength at the photosphere, since the field lines are 
supposed to open out with height even in the chromosphere 
(e.g. Giovanelh 1980). 



3.2 Boundary Conditions and Numerical Method 

Now we would like to explain the practical aspects of our 
method of constructing a unique coronal wind structure with 
respect to various input properties of the linearly polarized 
Alfven waves and the acoustic waves. In order to solve both 
the heating of the corona (energy transfer) and the accelera- 
tion of the solar wind (momentum transfer) consistently, our 



calculation is performed in a broad region from the bottom 
of the TR where temperature, Tp = % = Tin = 2 x 10*K, 
located at Hh = LOOSii© (e.g. Golub & Pasachoff 1997) to 
an arbitrary outer boundary at rout = 215-R0(= lAU). 

Only transonic solutions are allowed for the flow speed, 
V. Numerical integration of eas.CT. I20ll . I23II . 1251 . and 1261 1 
is carried out simultaneously from the critical point, r = ra, 
deflned below by Rosenbrock method, which adopts i mplicit 
spati al grids to keep the numerical stability (Press et alJ 
Il992h . When we perform the integration of the energy equa- 
tions 1251 1 and 1261 1. isothermal sound velocities for the pro- 



tons, a„ = 



and for the electrons, a. 



used. To start the integration, we set seven initial guesses. 
Per, flp^cr, fleer, {dap/dr)ci, {da^ / dr)cy: , ash,x,cr, and ash,||,cr 
(subscript, 'cr' denotes the 'critical' point), r^r and Hcr are 
automatically derived from conditions that the numerator 
and denominator of a transformed form of the momentum 
equation 1231 1. 

dt; _ a? Qsh,x«A(^^A + 3«) 1 
dr V 4i;x (t^A + v)v 



( 2 O^sh.X Vp^{vA + 3v) , 



4cx 



+ 



VA +V 
h,X«A 



^dA 
A dr 



da^ 
dr 



GMr. 



8rx,i(^ + i'A)2(l-ci/«i)2J 



(38) 



should become zero at 



where a 



I cr, VViiV^iV. U,p 

■ ■ 

{dv/dr)cr is also obtained by differentiating both the nu- 
merator and denominator of ea. l38t (De I'Hopital's rule; see 
e.g. Lamers & Cassinelli 1999). 

The integration is firstly carried out to outward direc- 
tion to satisfy following outer boundary conditions by tuning 
{dap/dr)cr and {da^ / dr)cY'- 



V 
V 



i^c,p(rout) = 0, 

i^c,c(7-out) = 0. 



(39) 
(40) 



Although the integration is to be carried on to r — + oo ideal- 
istically, we have to stop it at r = rout . The above conditions 
are requirements to keep further outward integration stably 
without divergent behaviors of T ^ or T ^ oo. When per- 
forming the outward integration, we have to be careful when 
solving eq. 1201 for the switch-on shocks, since it is only valid 
in the low-/3 plasma [va > Cs). (Note that the second term 
of the right-hand side includes 1/(1 — /v\).) The condition 
of VA > Cs holds in a range of r ~ (40 — 120)i?Q in our 
calculations. We firstly solve eg. 1201 . and once the condi- 
tion of the low-/3 breaks, we set ash,x = instead of solving 
eq. 12U1 . This ad hoc treatment can be justified, since the 
wave energy of the switch-on shocks is almost completely 
dissipated in the low-/3 region. 

After fulfilling the outer boundary conditions, eos. 1391 
and 1401 . the inward integration is followed. If our initial 
guesses are good, Tp and Tc approach each other as r de- 
creases. Therefore, to satisfy the condition of 

(41) 



we look for the correct Op ^r and flower, where nc(— 3 x 
lO^cm"^) is coupling density above which the thermal cou- 
pling between the protons and electrons is well- achieved. 
When n > ric is realized, the term for the energy ex- 
change between the protons and electrons dominates the 
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other terms in eqs. (1^ and The solar corona is treated 
as the one-component plasma in the inner part. In our calcu- 
lations, the two-component plasma is converted to the one- 
component plasma around 1.1 — 1.4iiQ. 

Then, we integrate ea. l|29|l inwardly instead of eas. l|25|l 
and 1261 1 to r = rin to satisfy the inner boundary conditions 
below: 

Fw,±(n„) = fw,±(r_L,d) = i^w,±,o (42) 
F„,||(n„) = i^w,||('-||,d) = i^w,||,o (43) 
T(n„) = Tin, (44) 

|Fc(ri„)|(~0) < jfcmaxl, (45) 

where -Fc.max in ea. H45|l is the maximum value of the down- 
ward conductive flux in the inner corona. The first and sec- 
ond conditions denote that the wave energy flux of the lin- 
early polarized Alfven waves (transverse mode) and acous- 
tic waves (longitudinal mode) must agree with the given 
values when the waves start to dissipate. r±,d and ry^d are 
locations where the waves start to dissipate. For the acous- 
tic waves, we set ry^d = ^© + 2 x 10^ km, focusing on the 
waves with ry ~ 100s excited at Rq + lO'^km in the corona 
by some dynamical events (S02). For the linearly polarized 
Alfven waves, we choose rx,d = ^in, which corresponds to 
the case that the waves are generated at the photosphere 
and start to dissipate from the bottom of the corona. If one 
takes the waves excited in the corona, larger rx,d should 
be adopted. For example, rx,d — 1-SRq for the sinusoidal 
waves with t± = 120s ( H2.2t . However, the dissipation of 

the switch-on shock trains in the low corona (r ~ 1.5Re) 
is not significant (Fig0J, hence, the choice of different rx,dS 
affects the results little. 

The third condition, ea. l44ll . is straightforward; the 
temperatures have to coincide with the fixed value at the 
inner boundary. The fourth condition, eg. 14511 . is the require- 
ment that the downward thermal conductive fiux should ra- 
diate away and become sufficiently small at the bottom of 
the TR (T = 2 X lO^'K), diminishing from its enormous value 
at the coronal base (T ~ 10® K). Practically, we continue 
calculations iteratively until Fc(rin)/-Fc,max < 1% is satis- 
fied. Note that thanks to this condition, the coronal base 
density, which is poorly determined from the obs ervations, 
does n ot have to be used as a boundary condition iHammerl 
ll982allbHwithbro(il98ll . The density at the coronal base or 
the TR is calculated as an output; a larger input energy in- 
creases downward conductive flux from the lower corona to 
the chromosphere, demanding a larger density in the coronal 
base and TR to enhance radiative cooling to balance with 
the increased conductive heating. 

Unless the above boundary conditions of eqs. 1421 - 14511 
are simultaneously satisfled, one returns to the outward inte- 
gration with preparing an improved set of initial guesses and 
iteratively determine the unique coronal wind structure on 
the given wave energy fiux. Physically, the conditions of the 
wave fiux (eqs. 1421 and 1431 ') regulates Ofsh.x.cr and Osh.ii.cD 
those of the temperature fegs. MH and 1441 ') regulates ttp ^.^ 
and fle.cr, and that of the conductive flux f eg. 1451 1 ) regulates 
Per. These relations guide the improvement of the respective 
initial guesses, though they are not independently approved. 
Finally, we can determine the unique coronal wind structure 



by iteratively improving the seven initial guesses to fulflll the 
seven conditions of eos. l39ll - 145ll . 



4 DEPENDENCY ON INPUT PARAMETERS 

We inspect dependency of the results on the input parame- 
ters, -Fw,x,o, F^,\\,o, T±, ry, /max and Bo to prepare for the 
next section in which we explore the model parameters ex- 
plaining the observed high- and low-speed solar wind. From 
now, we express period of the switch-on shock trains seen by 
the stationary observer just as t± instead of rx,i. We show 
dependencies of six characteristic properties, proton flux, 
(npii)iAu, at lAU, solar wind speed, wiau, at lAU, peak 
temperature of the protons, Tp^max, and electrons, Te.max, 
radiative flux, -Frad, integrated in the entire corona, and gas 
pressure, ptr, at location where T = lO^K in the TR. The 
physical interpretations of these six properties can be given 
as follows (e.g. Lamers & Cassinelli 1999): (npii)iAu is deter- 
mined by the momentum and heat input in the subcritical 
region, r < Vd. viau reflects the momentum and heat input 
in the supercritical region. Tp^max and Tc,max are determined 
by the heating from the wave dissipation and also by the 
coUisional coupling between the protons and electrons, ptr 
is essentially the same as the density at the inner bound- 
ary, frad reflects the density at the coronal base, since the 
radiative flux is oc for the optically thin plasma and the 
most of the radiation comes from the dense inner corona 
(< 1.5i?o). 

As for parameters, h± and fty, which control the energy 
distribution to the protons and electrons from the wave dis- 
sipation, we only present the results of the equipartition case 
adopting h± — /ly = 0.5. This must be reasonable for the 
N-waves, since they work only in the inner corona where the 
protons and electrons are sufficiently coupled. On the other 
hand, the approximation might be incorrect for the switch- 
on shock trains travelling in the outer region (r > 2Rq). 
More energy is injected to the much more massive protons, 
while the lighter electrons move across the shocks trans- 
parently without getting so much energy. Consequently, h± 
might be smaller, though it is inevitable to perform the ki- 
netic simulation to determine h± correctly. We have per- 
formed the calculation for several models by adopting the 
extreme value of h± = 0, corresponding to the situation 
that all the dissipated energy is supplied to the protons. 
Although these calculations yields larger proton tempera- 
ture and smaller electron temperature than the equiparti- 
tion cases, they give quite similar structures of the density 
and the wind velocity. 

Figure|5|illustrates dependency of the output properties 
on the input parameters in contours. The fiducial values are 
set as (Tw,||,o,-Fw,x,o,Tx,-ry,/max,-Bo) = (3.0 x 10^,2.0 x 
10^,120, 120,5,5), where Fw,x,o and fw.n.o in erg cm"^s~\ 
T± and Ty in second, and Bo in Gauss. Four panels on the top 
show dependence on 0.5 x 10* ^ F^,±,o, F^,\\,o ^ 5.5 x 10*, 
those on the middle show dependence on 60 ^ ''"-Li ''"II ^ 180, 
and those on the bottom show dependence on 1 ^ /max ^ 20 
and 1 ^ Bo ^ 20. The first column presents contours of 
(wpw)iAu, the second column presents those of viau, the 
third column shows those of Tp^max and Te.max, and the 
fourth column shows those of F^ad and ptr- Thick dots in 
panels on the top and bottom denote models for which we 
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Figure 2. Contours showing dependency of output properties of the solar corona and the solar wind on the input parameters. Panels 
on the top show dependency on -Fw,±,o S'liti ^w,||,Oi those on the middle show dependency on rx and ry, and those on the bottom 
show dependency on /max and Bq. Panels on the first column present (npj;)iAu(10*cm~^s~^), those on the second column present 
iiiAuC^m/s), those on the third column show rp,max(10®K) (solid lines) and Te,max(10®K) (filled contours in gray scale), and those on 
the fourth column show _Fi.ad(10^erg cm~'^s~^) (solid lines) and ptr(dyn cm~'^) (gray contours). Thick dots indicate the parameter sets 
for which we fail to integrate the energy equations due to too low temperature in the outer region. Dashed lines in panel in the second 
column on the bottom indicate constant -Bo //max =4, 2, & 1. 



cannot determine the unique coronal structure because we 
fail to integrate stably the second order derivative of the 
thermal conduction term due to too low temperature in the 
outer region. 

Four panels on the middle indicate that dependence 
on the wave periods is quite weak within a range of 60 ^ 
T±,T|| ^ 180. In particular almost no dependence results in 
with respect to ry because the most of the energy of the 
acoustic waves is lost as the downward thermal conduction, 
being unrelated to ry. Weak but finite dependences on t± 
are consistent with the results obtained by Nakariakov et 
al.(2000). 

Dependency of (npi')iAu can be understood in terms 
of the energy and momentum deposition in the subcritical 
region. Larger input of -FW.x.o and fw,||,o directly enhances 
the energy and momentum inputs to raise (npD)iAu (panel 
in the first column on the top) . Smaller Bo makes the switch- 
on shocks dissipate more quickly, since the shock strength, 
Q!ah,±, (eg. 1201 ). is scaled as oc oc -B,, (we neglect v 



in ea. (l34ll ) for fixed -Fw.±- As a result, small Bq increase the 
energy and momentum injection in the subcritical region to 
enhance (npi;)iAu (panel in the first column on the bottom). 

A panel in the second column on the top shows that 
higher «iau is achieved by models with smaller F^^y.o- This 
is because smaller input wave energy flux gives less dense 
corona and the smaller amount of coronal gas can be effec- 
tively accelerated to higher velocity. One of the most im- 
portant panels is that in the second column on the bottom 
which illustrates that vwv is well correlated with -Bo //max- 
Bo/ f max determines -By in the outer region of r > ri + cr 
where the flow tube expands almost radially (eg. 13711 ). In 
models adopting larger -Bo//max the acceleration and heat- 
ing occurs more in the supercritical region than in the sub- 
critical region because ash,± oc B^^^'^ , and larger «iau 
yields. Interestingly, the nice correlation betw een v\av and 
Bo/.f max is also obtained in recent observation iHirano et alJ 
l2003h . which might be explained by the similar discussion. 

Comparison of panels in the third column with those 
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in the first column show that correlation between To, max 
and (ripi')iAu is good. This indicates that T^max is mainly 
determined by the energy input in the subcritical region. On 
the other hand, Tp.max is regulated not only by the energy 
deposition but by the coupling with the electrons; too large 
input of energy in the subcritical region decreases Tp^max 
since more energy is transfered to the cooler electrons due 
to the larger density. 

Three panels in the fourth column illustrate that F^ad 
and pti are determined almost solely by F„^\\fi. The N- waves 
rapidly dissipate and the heating occurs mainly in the inner 
corona ( ~ 1.5Rq). The dissipated energy is mostly lost as 
the downward condu ction which finally radiates away in the 
TR llHammeilll982allbh . The radiative loss is determined by 
the density there according to oc p . Hence, the energy from 
the N-wave dissipation controls the density at the coronal 
base and the TR. In other words, larger input of the energy 
in the inner corona can heat the plasma deeply down to the 
chromosphere. 

Before closing this section, we summarize several im- 
portant points with respect to the dependency on the input 
parameters. 

(i) Jrad and pti, or equivalently the density at the coronal 
base, is mostly determined by -Fw,||,o- 

(ii) The solar wind speed, «iau, is well correlated with 
So//max, whereas it is anticorrelated with -^^,^1,0- 

(iii) Within a range of 60s ^ ''"-L)''"!! 180s the coronal 
wind structures are affected by t± and ry little. 

(iv) (npw)iAu, mass flux of the solar wind, depends on all 
the parameters. 



5 COMPARISON WITH THE OBSERVED 
HIGH/LOW-SPEED SOLAR WIND 

5.1 Structure of Corona and Solar Wind 

To test reliability of our process in the real solar corona, 
we study whether our model can account for the difference 
between the high-speed solar wind from the polar coronal 
holes and the low-speed solar wind in the equatorial region 
during the low solar activity phase. The aim of this section 
is not to show superiority of our process to other heating 
and acceleration sources but to explore possibilities whether 
our process can become one of the reliable mechanisms. Re- 
cently, two or three dimensional structures have been inves- 
tigated extensively by numerical simulations (e.g. Roussev 
et al.2003). However, the self-consistent treatment of the 
wave propagation is still difficult in such multi-dimensional 
calculations. Therefore, our one dimensional studies play a 
complementary role with them. 

As we saw in Fig|51 the variation of the wave periods, 
60s ^ TXjTii ^ 180s, affects the results little. Therefore, we 
look for the model parameters, -Fw.x.o, F^,\\,Oy Bo, and /max, 
for fixed t± = ry = 120s for simplicity. We summarize be- 
low how we determine the input parameters, Fw,x,o, ^^w,||,0) 
Bo, and /max, which reproduce the data best, based on the 
arguments in the previous section. 

(i) The density in the inner corona is regulated almost 
solely by Fw,||,o- The observed density in the inner corona 
( ~ 1.2i?0) can be explained by J'w.n.o — 3 — 5 x lO^erg 



cm~^s~^ for the less dense polar holes and F„,\\,o — 7 — 9 x 
lO^erg cm~^s~^ for the denser streamer region (Fig'OJ. 

(ii) Since the solar wind speed at lAU is well corre- 
lated with -Bo//max, larger -Bo(G)//max ~ 2 is re- 
quired for the high-speed wind {viAv ~ 700km/s) and 
smaller -Bo (G) //max — 1 is favorable for the low-speed wind 

{viAv ~ 400km/s). 

(iii) The observed velocity profile in the polar coronal 
holes shows high outflow speed in the low corona, while the 
low-speed wind appeared to be accelerated gradually (bot- 
tom panel of Fig|3J. Jumping to the conclusion, smaller Bo 
is required to reproduce larger v in the inner corona. Since 
the flow tube does not completely open to /max in the in- 
ner corona, smaller Bo results in smaller field strength, , 
there. Smaller anticipates larger ash,± for fixed -Fw.x, 
which let the transverse wave dissipate more in the inner 
corona. The input of more energy in the subcritical region 
enhances mass flux, pv. On the other hand, p in the inner 
corona is independently determined by -Fw.n.o- Consequently, 
a model adopting smaller Bo with identical Fw,±,o, F„,\\,Oj 
and Bo //max gives larger v in the inner corona. 

(iv) From the above restrictions of (i)-(iii), y 0, Bo, 
and /max are fixed. The last parameter, -Fw,x,o, is tuned to 
explain the density and the temperature in the intermediate 
region, 2Rq ~ r ~ 6Rq ; too much input of the transverse 
waves gives too high temperature and also too high density 
due to too large density scale height, and vice versa. 

We tabulate the adopted model parameters and the results 
of the physical properties of the corona and the solar wind 
in TabQ The structures of the corona and the solar wind in 
both models are presented in Figl^with the observational 
data. 

As shown in Tab^we have adopted larger -Fw,||.o by a 
factor of 20 in the low-speed wind model than in the high- 
speed wind model based on the condition (i). This directly 
leads to larger radiative flux in the low-speed wind model by 
a factor of 13, which is consistent with observational result 
that X-ray flux in the quiet-Sun in the low latitude exceed 
that in the polar coronal holes by more than an order of mag- 
nitude (e.g. §1 of Golub and Pasachoff 1999). Accordingly, 
the transverse waves become relatively important against 
the slow waves (-Fw,±,o/^V,||,o 3> 1) in the high-speed wind, 
while the longitudinal waves dominant (-F'w_x,o/-Pw,||,o < 1) 
in the dense low-speed wind. Table also shows that the 
comparable order of the mass flux is realized in the high- 
speed wind model in spite of smaller input of -Fw,±,o + ^'w,||,0) 
since smaller /max is adopted and the effect of the adiabatic 
expansion of the flow tube is reduced. 

From the restrictions (ii) and (iii). Bo and /max for the 
high-speed wind model is almost uniquely determined as 

(Bo(G),/max) = (2,1). We choose (Bo(G),/max) = (10,8) 

for the low-speed wind model which shows a little wider ac- 
ceptable ranges, 6 ~ Bo ~ 12 and 5 ~ /max ~ 10. It 
should be noted that Bo, the field strength at the bottom of 
the TR, is smaller than the photospheric value which can be 
observed, provided that the field lines open out with height 
in the chromosphere. Although this ambiguity remains, our 
results of Bo are qualitatively consistent with th e empirical 
determination bv lSittler fc Guhathakurtal (1 19991) giving sur- 
face field strength ~ 3G at the pole and ~ 15G at the equa- 
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Table 1. Models for the higli/low-speed solar wind; iv^x.Oi ^w,||.Oi ^^'^ ^rad in (10^ erg cm^^s^-'^), Bq in Gauss, (npt))iAU in 10^cm~^s~^, 
^'lAU in km/s, Tp^ 

max s-nd ^c,max in lO^K, and ptr in 10 ^dyn cm ^. Note that is field strength at the bottom of the TR and /max 
is also normalized there. 
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Figure 3. Comparison of our results of the high-speed solar wind 
(solid lines) and the low-speed solar wind (dashed lines) with the 
observed data. In top and middle panels, open symbols are data 
in the polar coronal holes during the low activity phase and filled 
symbols are data in the equator or mid-latitude region where 
the low-speed solar wind is supposed to originates from. Top : 
Distribution of electron density. Middle : Distribution of proton 
and electron temperature. Note that proton temperature is higher 
than electron temperature in both models, bottom : Outflow ve- 
locity of the solar wind. Shaded regio n is observational data in 
the streamer belt iSheele v et al."l997^. Data with error bars are 
outflow speed b ased on ob servation of the interplume region in 
the polar holes iTeriaca et al.ll200t}) 



tor. In a region of r ~ 2Rq the field strength (oc Bo/Zmax) 
in the high-speed wind is lager than that in the low-speed 
wind, which is opposed to the situation in the low corona, 
due to the flow tube expansion. 

The adopted value, /max = 1, for the high-speed wind 
model corresponds to the case of the radial expansion of 
the flow tube in the corona, which is consistent with the 
observed radi al structure around 1 .15 — 5.5Rq in the polar 
coronal holes (IWoo fc Habball999^ . However, we would like 
to point out that our /max is measured from the bottom of 
the TR, which may be smaller than the observed expansion 
factor normalized by the photospheric field strength. Our 
results coul d become consistent w ith a reported expansion 
factor ~ 4 llGoldstein et alJll99d) . if the flux divergence of 
a factor of 4 occurs mainly in the chromosphere. 

On the other hand, our analysis has given relatively 
large /max for the low-speed wind model, which means that 
the closed structures cover large fraction of the surface. 
There have been two origins considered with respect to the 
low-speed wind (see Wang et al.2000 for review). One is the 
acceleration due to the intermittent break-up of the cusp- 
shaped closed fields on the equatorial region (e.g. Endeve, 
Leer, & Holtzer 2003). The other is the acceleration in the 
open flux tube with large areal expansion in the low- and 
mid-latitude corona (e.g. Kojima et al.l999), which can be 
investigated by our modelling. Koiima ct al. ( 199!J) detected 
the low-speed wind with a single magnetic polarity, originat- 
ing from the open structure region located near the active 
region. They showed the rapid expansion of the flow tube, 
being qualitatively consistent with our result, whereas /max 
cannot be directly compared with the reported expansion 
factor ~ 50 due to the same reason described above for the 
high-speed wind model. 

The top panel of FigO shows electron density of 
the high- and low-speed wind model, overlayed with ob- 
servational data in the polar region (open circles are 

]); open stars 
ope n trian- 



SVMER/ SOHO data bv lWilhelm et al 



VVCS/ SOHO data bv lTeriaca et~ 



gles and sq uares are LA S CO / S OHO data bv lTeriaca et alJ 
(2003) and iLamv et alj il997t) ') and in the streamer re- 
gion (filled diamonds and triangles are CDS/ SOHO an d 
VYCS/ SOHO data at mid-latitude bv lParenti et al.l J200CD : 
fill ed pentagons are hASCO/ SOHO data at the equator 
bv iHaves. Vourlidas. fc HowardI ((20^)). As for the obser- 
vations in the polar holes, we present data in the interplume 
regions whe n the plume and int erplume are distinguished, 
according as lTeriaca et alJ (l2003l) reported a strong evidence 
in favor of the interplume regions as sources of the high- 
speed wind. The panel exhibits that our results can explain 
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the trends of the observed electron density well in both mod- 
els. 

The middle panel of Figgl compares the results of 
proton and electron temperature with the ob served pro- 
ton te mperature (data with error bars are by lEsser etlil] 
^1999^: open squares are VYCS / SO HO observation by 
lAntonucci. Dodero. fc Giordanol i2000f) l and electron tem- 
perature (open circles are C DS/ SOHO observation by 
iFludra. Del Zanna. fc Bromaed 1(1999)) in the polar holes 
and the observed electron temperature in the mid- 
latitude streamer (fil led diamonds and triangles are 
CDS,\JYCS/ SOHO by iParenti eTai] j200Ci) ). Please note 
that the proton temperature is higher than the electron tem- 
perature in both models. The results of the proton and elec- 
tron temperature for the high-speed solar wind show rea- 
sonable fits to the data. Our results appeared to be slightly 
lower than the observed kinetic temperature of the protons 
by Antonucci et al.(2000; open squares). These data indi- 
cate temperature perpendicular to the radial magnetic field. 
Although our model consider the isotropic distribution of 
the temperature, it is expected that the heating by the fast- 
shocks is anisotropic and more effective in t he perpendicular 
direction to the field line jLee fc Wull200(]|) . Therefore, our 
model could give better fitting to the data if taking into 
account the anisotropic temperature distribution. On the 
other hand, the temperature for the low-speed wind can- 
not reproduce the observed high temperature in a region of 
1.3Rq ~ r ~ 1.6i?0 . It is difficult to reproduce the obser- 
vation only by adopting larger t± and ry which give slower 
dissipation. Therefore, we have to say that some other heat- 
ing mechanisms are required in this region. 

The bottom panel of Fig|3| presents outffow velocity of 
the solar wind. The result of the outffow velocity for the 
low-speed wind is well inside of the obtained range by the 
observation of about 65 moving objects in the streamer. On 
the other hand, our result for the high-speed solar wind is 
marginally consisten t with the data of the lower edge by 
iTeriaca et alJ ll2003h . A cooperation with an extra but less 
dominant acceleration process would explain the data com- 
pletely. 



5.2 Propagation of Shock Trains 

Figure 0] compares the N-waves and the switch-on shock 
trains in the high- and low-speed wind models. It shows 
that the dissipation of the N-waves is much more rapid than 
that of the switch-on shock trains in both models. Qsh,|| and 
F^,\\ is drastic decreasing functions of r, while Ofsh.x is more 
or less constant and -FW.x decreases gradually. Consequently, 
dissipation length, |F,„^||/V • F^,\\\, of the N-waves is very 
small, ~ O.OI-R0 in r ~ 1.5Rq. These results indicate that 
the dissipation of the slow MHD waves is a powerful heating 
source in the inner corona, which is clearly illustrated in the 
top panel of Fig|^ comparing heating per unit mass, Q±/p 
and Qii /p. The role of the dissipation of the slow MHD waves 
is important even in the high-speed wind model, adopting 
the relatively small -Fwji.O) to sustain the moderate density 
at the coronal base through the downward conduction. The 
top panel of Fig^Jalso shows that variations of Qsh,|| are re- 
semble each other in spite of the large difference of the input 
energy ffux by a factor of 20. This is because the difference 



Shock Trains in H ig h/Low — S peed Wind 
1 2 5 10 20 

1 r . . . ^ 




Figure 4. Comparison between the switch-on shock trains and 
the N-waves in the high-speed wind (solid lines) and low-speed 
wind (dashed lines), top: Variation of shock amplitudes, cish,± 
and Osh,||i against radius, t/Rq. middle: Variation of wave energy 
flux, -F„^x a-nd ^w,||- bottom: Variation of dissipation length. To- 

tal dissipation length, |-= — ^^-^ (thick lines) is also 

V • (i'm,± + -Tw,!! ) 
displayed with dissipation length for the switch-on shock trains, 

V' — tJi^ — I, and for the N-waves, — — |. 
* w,_L V ■ ■'^w,! 

of the wave energy flux, _F„ y ~ pct^h ||'^s, mainly attributes 
to variation of p; larger (smaller) input of ^Iv.n.o requires 
larger (smaller) density through the energy balance between 
the downward thermal conduction and the radiative loss as 
discussed so far. 

The heating from the N-waves is soon dominated by the 
heating from the switch-on shocks at ~ I.IRq in the high- 
speed wind model and 1.2_Rq in the low-speed wind model. 
The heating from the switch-on shock trains keeps till a 
more distant region, which is seen in the top panel of Fig|^ 
Figure |1] further shows that the dissipation of the switch- 
on shock trains is more rapid in the low-speed wind model, 
since the shock amplitude, Qsh,± , is larger in the inner corona 

(the top panel). This is because ash,±(oc p^/^B^^ F^^±_) 

becomes systematically larger in the dense low-speed wind 
with smaller -Bo//max. 

The dissipation of the switch-on shock trains also con- 
tribute directly to the acceleration by the gradient of the 
wave pressure (eg. 1301 ). The bottom panel of Fig|K|presents 
the acceleration per unit mass, dpv,^±/ pdr, from the switch- 
on shocks compared with the acceleration from the gas 
pressure gradient, dp/pdr. Note that the acceleration from 
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Figure 5. Comparison of heating and acceleration in the high- 
speed wind (Solid lines) and low-speed wind (dashed lines), fop: 
Heating per unit mass from the switch-on shock trains, Qx/p, 
and from the N-waves, Q\\l p. bottom: Momentum input per unit 
mass from the switch-on shock trains, dp„ j_/{pdr), is compared 
with that from the gas pressure, dp/{pdr). 



the N-waves is assumed to be zero as it is negligible. The 
panel indicates that the acceleration by the switch-on shock 
trains dominates that by the thermal pressure in a region of 
3 — lOi?0 in the high-speed wind model, while the gas pres- 
sure dominates in the entire region in the low-speed wind 
model. In summarize, the switch-on shock trains contribute 
both to the heating and to the acceleration in the high-speed 
wind, while the low- wind wind is basically driven by the 
thermal pressure of the coronal plasma which is originally 
heated by the dissipations of the switch-on shock trains and 
the N-waves. 

The results shown in the bottom panel of Fig|l] allows 
us to check the validity of the conventional assumption of 
the constant dissipation length. Total dissipation lengths, 
|(Fw,x + F„,||)/(V ■ (Fw,_L + -Fw,||))|, (thick lines) are at 
first subject to those of the N-waves in the low corona and 
eventually to those of the switch-on shock trains in the upper 
region ( ~ 1.3Rq). Because of the transition, the total 
dissipation lengths rapidly increase. Even in the outer region 
they are not constant; they vary 2 — 6Rq in the high-speed 
wind model and 0.4 — 2Rq in the low-speed wind model. 
Therefore, the assumption of the constant dissipation length 
is not adequate for our processes. 

We would like to discuss the amplitudes, ash,± and 
aah,||, in connection with the observed non-thermal broad- 
ening. Velocity amplitude of the transverse mode in the 
high-speed wind model becomes < 5v± >= -^Qsh,±i'A ~ 

1000(km/s) X 0.15/^3 ~ 90(km/s) at its maximum 
level, where a factor, -4j, comes in for the trains il- 
lustrated in Fig0 On the other hand, the off-limb ob- 



servation in the polar coronal hol es gives non-thermal 
broadening of Svoba — 30 — 50k m/s iBaneriee et al.l[l998l : 
iDovle. Teriaca. fc Baneried 119991) . which shows that our < 
Sv± > is not consistent with Svoha- However, since uncertain- 
ties of the projection effects might affect Svoha, it would be 
hasty to conclude that our model is ruled out by the observa- 
tions. Velocity amplitude of the longitudinal mode is as large 
as (5ii|| = -^Q!sh,||Cs ~ 50km/s at the coronal base in both 
models, and it decreases to (5ii|| < lOkm/s ai r — 1.3Rq. 
Nonther mal broadening of 2 0-40km/s is observed in the so- 
lar disk l)Erdelvi et al.ll9'9^ . Our results are consistent with 
the observed result except at the coronal base. It is also un- 
certain whether the discrepancy at the coronal base is real 
or not, since the same ambiguity comes into the observations 
of the longitudinal mode. 



6 SUMMARY 

We have investigated the coronal heating and the acceler- 
ation of the high- and low-speed solar wind in the open- 
field regions by the dissipation of the shock trains. The 
small reconnection events as well as the surface granula- 
tions excite MHD waves not only at the photospheric level 
but in the corona. Among these waves, the linearly polar- 
ized Alfven waves and acoustic waves steepen through the 
upward propagation and travel as the switch-on shock trains 
and N-waves respectively. 

Firstly, we derived evolutionary equations for the am- 
plitudes of the shock trains under the WKB approximation 
by assuming that the dissipation occurs only through the 
shocks. While the assumptions are supposed to be reason- 
able for the longitudinal mode, we have to be careful in 
dealing with the transverse mode because other types of the 
dissipation mechanisms may be important as well as the 
WKB approximation becomes marginal in the low corona 
for the waves with periods of a few minutes. In the fu- 
ture studies, it is require d to treat the transvers e wave s 
dynamica lly as taken i n lOfm an fc Davila I il997l Il998l) : 
GraDoin, Leorat. fc HabbaD(l2002l) . whereas they consider 
the wave propagation under the one-fluid approximation. 

We determined coronal wind structure consisting of the 
protons and electrons from the bottom of the TR to lAU 
by solving the wave equations simultaneously with the mo- 
mentum and energy equations. Studies on the dependency 
on the input parameters give the two important results: 

(i) vau has nice correlation with the field strength in the 
outer region which is proportional to Bo //max. 

(ii) The density at the coronal base is almost solely de- 
termined by -Fw.ii.o- 

In the low-speed wind, the acoustic waves play a role to 
keep the dense corona and smaller Bo //max — 1 is required 
to hold the wind speed slow. On the other hand, in the 
high-speed wind the transverse waves become relatively im- 
portant and larger -Bo//max ~ 2 is favorable. The early on- 
set of the acceleration in the high-speed wind also requires 
the weak field strength in the low corona. This accordingly 
leads to the radial expansion of the flow tube in the corona, 
though it is not inconsistent with the observed expansion 
factor ~ 4 in the polar holes if the flux divergence occurs in 
the chromosphere. 
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Our model has reproduced the overah trend of the high- 
and low-speed solar wind except the observed high temper- 
ature in the mid-latitude corona connected with the low- 
speed solar wind. Therefore, our processes could be reliable 
heating and acceleration mechanisms, whereas contribution 
from other heating sources are required in the low-speed 
wind region. 
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